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ABSTRACT 

We present new (April 2000) MERLIN observations of the H 2 0 masers located near the protostar 
Cepheus A HW2. The MERLIN observations detect many of the structures found in earlier (1996) 
VLBA observations of Torrelles and collaborators, and the changed positions of these structures are 
compatible with the VLBA proper motions and astrometric uncertainties. The radius of curvature of 
the R4 structure of maser arcs appears to have grown by a factor of two, and the displacement of the 
arcs between 1996 and 2000 are compatible with expansion about a common center. In addition, the 
MERLIN observations detect red-shifted masers not previously found; taken with the newly discovered 
masers, the R4 structure now resembles patchy emission from an elliptical ring. We demonstrate that a 
simple bow-shock model cannot simultaneously account for the shape and the velocity gradient of the 
R4 structure. A model involving a slow, hydromagnetic shock propagating into a rotating, circumstellar 
disk better describes the maser spot kinematics and luminosities. In this model, the central mass is 3 Mg, 
and we demonstrate that the mass of the disk is negligible in comparison. The expansion velocity of the 
post-shock gas, ~ 5 km s , is slow compared to the shock velocity, vs ~ 13 km s -1 , suggesting that the 
post-shock gas is magnetically supported with a characteristic field strength of ~ 30 mG. We speculate 
that the expanding maser rings R4 and R5 may be generated by periodic, instability-driven winds from 
young stars that periodically send spherical shocks into the surrounding circumstellar material. 

Subject headings: ISM: individual (Cepheus A) — stars: formation — stars: planetary systems: 
protoplanetary disks - - stars: pre-main sequence — masers 


1. INTRODUCTION 

Interstellar H 2 0 masers trace warm, dense molecular 
gas associated with young stellar objects (YSOs) and star¬ 
forming regions (e.g., Genzel & Downes 1977; more recent 
examples include Claussen et al. 1998; Furuya et al. 2000; 
Rodriguez et al. 2002). In contrast with thermal tracers, 
H 2 0 masers are bright, compact radio sources, and can 
be detected on very long baselines, corresponding to sub- 
milliarcsecond (mas) resolution. Using aperture synthesis 
techniques, one can map both the geometry of the maser 
spots and trace their very fine proper motions. Adding 
radial velocity information derived from the Doppler shift 
of the maser line results in a detailed picture of the maser 
geometry and kinematics as a measure of the dynamics 
very close to the YSO (e.g., Genzel et al. 1981a; Genzel 
et al. 1981b; Schneps et al. 1981; Gwinn, Moran, & Reid 
1992). 

Most of the known interstellar masers trace gas lo¬ 
cated near the base of molecular outflows on scales of 
10s 100s of AU (for reviews, see Reid & Moran 1988; 

Elitzur 1992) The association of masers with outflows has 
been unambiguous. The maser spots (i.e., individual, 


unresolved maser sources) align with molecular outflows 
mapped in millimeter-wave emission lines, or are associ¬ 
ated with Herbig-Haro objects. The kinematics of the such 
masers also follow the sense of the motion of the outflow 
on larger scales. The proper motions are usually small 
compared to the molecular outflow velocities, suggesting 
that H 2 0 masers trace the expanding shock front between 
the outflow and the ambient ISM (Claussen 2001). 

There are, however, a handful of studied interstellar 
H 2 0 maser sources that appear to be associated with cir¬ 
cumstellar disks or rings rather than molecular outflows 
(Matveenko 1987; Cesaroni 1990; Fiebig et al. 1996; Tor¬ 
relles et al. 1996; Berulis, Lekht, & Mendoza-Torres 1998; 
Torrelles et al. 1998a; Torrelles et al. 1998b; Hunter 
et al. 1999; Shepherd & Kurtz 1999; Matveenko & Dia¬ 
mond 2000; Patel et al. 2000; Lekht & Sorochenko 2001). 
These disks are thought to be the remnants of the molec¬ 
ular cloud cores out of which the stars formed (e.g., Shu, 
Adams, & Lizano 1987). The nature of these disks ulti¬ 
mately bears on the origin of planets, which are generally 
thought to form out of condensations in the circumstellar 
disk (e.g., as reviewed in McCaughrean 1997). 

Recently, Torrelles and collaborators (Torrelles, et al. 
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1996; (Torrelles, et al. 1998a) reported evidence for H 2 0 
masers in a protoplanetary disk around the continuum ra¬ 
dio source Ceph A HW2 (Hughes & Wouterloot 1984a). 
The radio continuum source resolves into a jet of ther¬ 
mal free-free emission oriented along PA 44° (Rodriguez 
et al. 1994). Torrelles et al. used the VLA to map the 
maser emission and found a compact distribution of spots, 
roughly 300 AU (0"4) in extent (assuming a distance of 
725 pc: Johnson 1957). The maser spots coarsely align 
at right angles to the jet axis, i.e., are better aligned with 
the predicted position angle of a protoplanetary disk than 
with the thermal jet. In addition, there is a radial veloc¬ 
ity gradient across the distribution of maser spots that is 
compatible with disk rotation. 

Their subsequent 1996 VLBA observations resolved the 
apparent disk into a more complex, filamentary arrange¬ 
ment of maser spots. The maser spot filaments (called 
R1 R5) show seemingly independent proper motions and 
may well be associated with unseen protostars in the neigh¬ 
borhood of HW2 instead of a large-scale disk associated 
solely with HW2 (Torrelles et al. 2001a; Torrelles et al. 
2001b). In this work, we present new MERLIN observa¬ 
tions of the H 2 0 masers located near HW2. Our main 
result is the discovery of an expanding ring of maser spots 
associated with region R4 (using the naming convention 
of Torrelles et al.). There is a velocity gradient along the 
major axis, which is most simply explained by rotation. 
We consider the possibility that the R4 masers arise from 
a gaseous disk surrounding a forming star. 

This paper is organized as follows. §2 describes both 
the new MERLIN and archival VLA observations of the 
HW2 H 2 0 masers, and §3 presents the results of these ob¬ 
servations. §4 describes two shock models for the HW2 
R4-A maser arcs and argues in favor of a model involving 
a shock wave propagating into a rotating, gaseous disk. 
We also present some constraints on the properties of the 
pre-shock gas based on the shock models. §5 considers the 
proper motions of the other arcuate maser structures, R1 
R5 with the goal of placing a limit on the proper motion 
of the rotation center of the R4 masers. §6 summarizes 
the results of the observations and shock modeling and 
also briefly discusses the possible origin of the R4 masers. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. MERLIN Maser Data 

On 9 April 2000, we used the 5-element Multi-Element 
Radio-Linked Interferometer Network (MERLIN 1 ) tele¬ 
scope, based at the Jodrell Bank Observatory, to observe 
the Ceph A HW2 region (pointing center: a[1950] = 22^ 
54 m 19.0592 s , 5[1950] = 61°45'46.51"). The central fre¬ 
quency was tuned to the 6i6 - 5 2 3 H 2 0 maser transition 
(rest frequency: 22235.079 MHz), offset to the systemic 
velocity of the Ceph A region, Vlsr = —12 km s -1 . The 
observations were made in spectral line mode with 256 
channels and a channel spacing of 0.21 km s -1 . The total 
on-source integration time was 11.5 hrs. 

Calibration involved initial editing and bandpass cali¬ 
bration using the “dprocs” routines developed at Jodrell 
Bank. We then used AIPS to compute an improved band¬ 
pass calibration, as well as atmospheric and instrumental 


phase and amplitude corrections against the position of 
the calibrator point source J2302+640, whose coordinates 
are known to a precision of about 12 mas (Patnaik et al. 
1992). The integrated maser spectrum derived from the 
calibrated visibility data is plotted in Figure 1. 

After calibration, the data were transformed into images 
of the sky using the AIPS task IMAGR, which performs 
both the requisite Fourier Transform and a “CLEAN” de- 
convolution (Hogborn 1974; Clark 1980). To achieve the 
optimum angular resolution, the visibility data grid was 
given uniform weight during the Fourier Transform. The 
resulting angular resolution (“clean beam”) was 8 mas cir¬ 
cular. The RMS noise in signal-free channels is 25 mJy 
beam -1 , comparable to the expected thermal limit of ~ 20 
mJy beam -1 . However, channels containing bright maser 
emission were always dynamic range limited at a ratio of 
typically 400:1 (signal: noise). 

Each channel image was then examined for maser signals 
of signal-to-noise ratio (SNR) > 5, from which we compiled 
a database of maser positions and velocities. The positions 
of the MERLIN maser spots are plotted in Figures 2 and 
3. For convenience, and unless otherwise specified, offset 
coordinates are referenced to the position of the bright¬ 
est spot associated with the R4 maser arc. The reference 
position is a(J2000) = 22 h 56 m 17^97792, Dec (J2000) = 
+62° 01'49 , . / 4421. 

2.2. VLA Continuum Data 

We also obtained archival VLA continuum observations 
of Ceph A HW2 to compare with the new MERLIN H 2 0 
maser data. The data were originally obtained and pub¬ 
lished by Torrelles et al. (1998a). We re-reduced the data 
following standard techniques within the AIPS software 
package. The flux calibration and astrometry of the re¬ 
sulting image compare very well with those presented by 
Torrelles et al. 

We converted the radio continuum image to J2000 co¬ 
ordinates using the AIPS task “REGRD.” In performing 
the coordinate precession, we took care to correct for the 
fact that the B1950 VLA phase calibrator positions are 
given for equinox 1979.9 (e.g., Muxlow et al. 1995). The 
continuum image is included as an overlay in Figure 2. 

3. RESULTS 

Figure 1 shows the spectrum of the calibrated visibility 
data. The spectrum is compatible with previous observa¬ 
tions in the broadest sense: there are many narrow, bright 
maser features spread from Vlsr ~ — 30 km s -1 to ~ 0 km 
s -1 . The spectrum does not agree in detail, however, ow¬ 
ing to the rapid variability of the masers associated with 
Ceph A; factor-of-two variations occur on timescales as 
short as a few days (Rowland & Cohen 1986). 

The distribution of maser spots relative to the posi¬ 
tion of the continuum source HW2 is plotted in Figure 2, 
and the distribution of maser spots in both position and 
(grayscale-coded) velocity is displayed in Figure 3. Ex¬ 
cluding the R5 maser group, the spots spread over ~ 1" 
east-west and ~ 0.5" north-south. Looking at the data 
qualitatively, there is a broad trend of increasing veloc¬ 
ity (blue-shifted to red-shifted) from east to west. These 
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Fig. 1.— Plots of the MERLIN total power spectrum of the H 2 O masers in the region of Ceph A HW2 and HW3. This spectrum was 
derived from the calibrated visibility ( u , v) data rather than the channel images. The top panel is the spectrum scaled to show the brightest 
emission (which is associated with HW3), and the bottom panel is the same spectrum scaled to emphasize the fainter masers. 
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Fig. 2.— Comparison of the MERLIN maser positions (small crosses) with the 22 GHz VLA radio continuum image (contours). The 
contour levels are 0.7, 1.2, 2.1, 3.5, and 6.0 mJy beam -1 . The VLA beam size, plotted in the lower left corner, is 0.08 // . The location of two 
newly discovered, faint continuum sources, VLA-mm and VLA-R5 (Curiel et al. 2002), are plotted as eight-point stars. 
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Fig. 3.— Plot of the HW2 maser positions and velocities based on the new MERLIN data. The axes are sky offset positions relative to the 
brightest maser within the MERLIN data. Each spot is shaded according to the radial velocity scale given in the bottom panel. Left panel: 
large-scale field showing all of the MERLIN detections over the HW2 region. Right panel: close-up field showing the details of the R4 region. 


results compare favorably with the previous VLA observa- 
tions of H 2 0 masers in HW2 (Torrelles, et al. 1996) and 
the arcsecond-scale velocity gradient observed in SiO (2 —> 
1) (Gomez et al. 1999). The improved spatial and veloc¬ 
ity resolution of the MERLIN data, however, reveals finer 
structure than could be mapped by the VLA. MERLIN 
resolves several arcuate groups of masers, including the 
R2 - R5 groups originally revealed by VLBA observations 
(Torrelles, et al. 2001a; torrelles, et al. 2001b). 

Torrelles et al. (2001a) presented details of only the 
blue-shifted masers associated with R4, and they measured 
a small proper motion away from the center of curvature. 
Accordingly, they interpreted the U-shaped structure as 
a bow shock caused by an outflow originating from a star 
near (in projection) to HW2 proper. The conditions within 
shocked molecular gas favor the generation of H 2 0 masers 
(Elitzur, Hollenbach, & McKee 1989; Kaufman & Neufeld 
1996), lending some support to this scenario. The MER¬ 
LIN observations reveal a structure that is nearly ring- 
shaped (Figure 3). The velocity gradient of these brighter 
masers runs counter to the larger-scale gradient of the 
HW2 region: velocities increase from the northwest to the 
southeast along PA 140° (see the fitting analysis below and 
in Table 2). This velocity gradient and elliptical geome¬ 
try suggest an alternative explanation: perhaps these R4 
masers trace a rotating ring of molecular gas surrounding 
the unseen protostar. The ring is probably only an annu¬ 
lar segment of a more extensive circumstellar disk. In the 
following section, we evaluate these two models in turn. 

4. MODELS FOR THE R4 MASERS 


4.1. Model 1: Expansion into a Static Medium (Bow 

Shock) 

The R4-A masers, those making up the northwestern 
arc of R4 (see Figure 3), present the most clearly defined 
structure of the region. A bow shock would naturally ex¬ 
plain the U-shape of the arcs (e.g., Raga & Bohm 1985; 
Raga 1986; Furuya et al. 2000) and the expansion of the 
R4-A structure, but it is not clear that bow shocks should 
produce the observed velocity gradient along the appar¬ 
ent major (outflow?) axis. The U-shape of a bow shock 
results from edge-brightening, which for masers would re¬ 
quire radial velocity coherence for effective amplification. 
The challenge for the bow shock model, then, is to have the 
resulting line-of-sight velocity coherence naturally provide 
both the arcuate shape of the maser spot distribution and 
the observed radial velocity gradient along the symmetry 
(outflow) axis (see Figure 3). 

We created a simple, parabolic shock front model to test 
this scenario. Following the formalism of Hartigan, Ray¬ 
mond, & Hartmann (1987) and Hartigan, Raymond, & 
Meaburn (1990), the geometry of the shock front is given 
by: 

z = a(x 2 + y 2 ) (1) 

where a is a shape parameter, z is the position coordinate 
along the outflow axis, and x and y are the position co¬ 
ordinates orthogonal to the outflow axis; the sense of the 
coordinates is illustrated in Figure 4. As gas clouds fall 
into the shock, they are accelerated only along the direc¬ 
tion normal to the surface of the shock. We will assume for 
this model that any bulk motions of the gas clouds are neg- 
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Fig. 4.— A sketch of the coordinate conventions used to analyze the bow shock model described in §4.1. 


ligible compared to their motion into the shock front. To 
evaluate the kinematics of the maser spots, which trace the 
post-shock gas, we need first to decompose the pre-shock 
velocity into components perpendicular to the shock front 
and parallel to the shock front. The post-shock velocities 
can then be estimated using the Rankine-Hugoniot condi¬ 
tions. A more rigorous treatment would take into account 
hydromagnetic effects. At this stage, however, the goal is 
to evaluate coherence effects, and a more complex model 
is not warranted. 

For the purpose of illustration, we assume isothermal 
(radiative) conditions for the shock. The actual conditions 
are probably closer to isothermal than adiabatic, because, 
in order to produce H 2 O masers, the pre-shock gas must 
already be dense, n (H 2 ) > 10 7 cm -3 (Elitzur, et al. 1989; 
Kaufman & Neufeld 1996). Based on the proper motion 
of the maser arc (see the discussion in § 4.2), the shock 
speed must be at least 13 km s , and so the immediate 
(adiabatic) post-shock temperatures are > 10 3 K (McKee 
& Hollenbach 1980). Under these conditions, the cool¬ 
ing time is only a few days, much shorter than the age of 
the maser arc. For the sake of completeness, it should be 
noted that an adiabatic shock model produces the same 
conclusions regarding radial velocity coherence; all that 
changes between the isothermal and adiabatic models is 
subtle, morphological details of the apparent structure of 
the shell. 

In the reference frame of the observer (i.e., at rest with 
respect to the stationary, pre-shock gas), the post-shock 
velocity is 

V = u'j_ + u'|| - Vs Z 

° l + 4aUP [2 ° ( ** + »y ) -z| (2) 

where u' ± and uj| are the velocity components of the post¬ 
shock gas in a coordinate system defined by the shock front 
(see Figure 4), vs is the shock speed, a is the shape pa¬ 
rameter defined above, and R = \/ x 1 + y' 2 . Thanks to the 
symmetry of the problem, we can take the observer to be 
located at some distance away from the shock front but 


confined to the ( x , z) plane (positioning in y is equivalent 
to an arbitrary rotation in PA). If the sight-line makes an 
angle i with outflow axis such that i = 0 is a pole-on view, 
the radial velocity of the post-shock gas is 

Vs 

VR=v o + -- 9 „ (2ax sin i - cos i) (3) 

1 + 4 a z R z 

where vq is the systemic velocity of the outflow source. 
For the purposes of simulating the maser arc, we esti¬ 
mated the velocity and shape parameters of the shock as 
follows. The maximum velocity of the R4 maser spots 
is ~ 7 km s -1 relative to systemic. Again assuming a 
fully radiative shock, deprojection gives the shock veloc¬ 
ity: vs « 7sec* km s _1 (in fact the proper motion argues 
for a shock speed of about 13 km s , however, the magni¬ 
tude of shock speed does not affect the conclusions). The 
projected shape of the shock on the sky (a 1 ) is related to 
the intrinsic shape ( a ) according to a = a' esc i (Harti- 
gan, et al. 1990). Based on a least-squares fit to the R4-A 
maser spot positions, we find a ~ 44 csc*(mas _1 ), which 
we used in our numerical models. 

We applied Equation 3 and the constraints on Vs and a 
to a set of models in an effort to mimic the properties of 
the R4-A maser arc; the results are shown in Figure 5. In 
these models, the only free parameter is the inclination of 
the outflow axis with respect to the line of sight. We as¬ 
sumed that masers would arise only within the Mach angle 
(i.e., where the perpendicular speed of the shock front is 
supersonic). We further assumed that the masers occupy 
a thin shell, appropriate for the short cooling time. The 
masers are preferentially amplified when the velocity dif¬ 
ference between the near and far sides of the thin maser 
shell is small. 

The model maser regions plotted in Figure 5 are appro¬ 
priate for a maximum velocity difference of 2 km s _1 , cho¬ 
sen to provide a reasonable match to the spread of maser 
spots around the R4-A arc. Relaxing either the thin shell 
assumption or the velocity difference criterion broadens 
the projected shape of the maser region, but the resulting 
maser spot velocity gradient along the outflow axis would 
be unaffected. From inspection of Figure 5, it is clear that 
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Fig. 5. — Results of the bow shock model. From top to bottom, each panel displays the bow shock model for inclinations i — 30°, 45°, & 
60°, respectively. Grayscale contours represent the radial velocity of the near side of the bow shock in steps of 1 km s —1 from systemic (light 
gray) to —8 km s —1 (black). The heavy, dark contour traces the region of the shock where the velocity difference between the near-side and 
far-side is less than 2 km s —1 and where the Mach number exceeds unity. Edge-brightening by velocity coherence would produce a maser 
distribution similar in shape and radial velocity as the region filling the heavy contours. 
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coherence effects naturally lead to edge-brightening and 
therefore can explain the U-shape of the R4 maser dis¬ 
tribution. The i = 30° model reasonably reproduces the 
extent of the R4-A arc. On the other hand, this simple 
bow shock model cannot reproduce the radial velocities 
of the maser spots. It seems that edge-brightening pref¬ 
erentially amplifies regions with radial velocities near the 
systemic velocity. In contrast, the maser spots show a 
systematic gradient of ~ 7 km s -1 along the inferred out¬ 
flow axis with the maximum velocities near the apex of 
the U-shaped distribution. The kinematics of the maser 
spots are not well-described by this simple model of a bow 
shock propagating into a static medium, and it seems clear 
that the pre-shock gas must be undergoing some system¬ 
atic bulk motion to explain the velocity gradient of the 
post-shock gas (i.e., the masers). 

It is worth pointing out that changing the details of the 
shock model do not significantly affect the results. In any 
model involving a bow shock, the post-shock gas expands 
along and away from the outflow axis. Sight-lines nearer 
the outflow axis intercept gas approaching the observer on 
the near side and gas receding from the observer on the 
far side. Velocity coherence is therefore poor near the out¬ 
flow axis. Sight-lines nearer the projected edge of the bow 
shock intercept gas moving more nearly in the plane of 
the sky whether on the near or far side. Radial velocities 
are therefore more coherent nearer the projected edges of 
a bow shock and lie close to the systemic velocity. 

4.2. Model 2: Expansion into a Rotating Medium (Disk) 

This model expands on the previous one and assumes 
that the pre-shock gas rotates around the source of an 
outflow or blastwave, presumably a protostar or a newly 
formed star. Whether the outflow is collimated or spher¬ 
ical, the rotational motion would be tangential to the 
strongest part of the shock front. The post-shock gas 
therefore largely retains the rotational component of mo¬ 
tion and picks up an additional expansion component. 
There remains, however, the issue of the arcuate struc¬ 
ture. Coherence effects only come into play if we observe 
the gas very near the equatorial plane of rotation, but an 
arbitrary geometry would not produce arcuate structure 
in this case. It seems likely, then, that the masers arise 
from dense, post-shock chimps of gas, without regard to 
larger-scale velocity coherence. 

The simplest geometry that can produce the arcs is a 
spherical shockwave expanding into an inclined, rotating 
disk. (Supposing an asymmetric shockwave adds orienta¬ 
tion parameters that we would be unable to constrain.) 
The blastwave shock-heats molecular material in the disk, 
leaving behind a temporary, patchy ring of maser emission 
(Elitzur, et al. 1989). To estimate the properties of the 
maser geometry subject to this blastwave and disk inter¬ 
pretation, we applied a tilted-ring model akin to that used 
for HI galactic velocity fields (Begeman 1989). To simplify 
the model and reduce the number of free parameters, we 
made the following constraints or assumptions. (1) The 
mass of the disk is much less than the mass of the central 
protostar, and therefore the rotation curve is Keplerian. 
(2) The position angle and inclination of the disk does not 
change significantly with radius and time (that is, there is 
little warping or precession between the 1996 and 2000). 


(3) The ring expansion velocities may have changed be¬ 
tween epochs, perhaps as a result of mass-loading as disk 
material is swept up. (4) The central protostar has not 
significantly accelerated between epochs such that the sys¬ 
temic velocity of the disk has not changed. To accommo¬ 
date these assumptions in practice, the model comprises 
four inclined rings, one for each epoch (three for the 1996 
VLB A observations). Each ring is constrained to have 
the same inclination and position angle, a common cen¬ 
tral mass to determine the rotation speed of each ring, 
and two expansion velocities, one for the 1996 epochs and 
the other for the 2000 epoch. To be clear, the expansion 
velocities are based on a purely kinematical fit to the radial 
velocities of the masers. In the context of this model, the 
expansion velocities measure post-shock velocities in the 
disk gas rather than the motion of the blastwave proper. 

Because the VLBA astrometry is referenced to a differ¬ 
ent sky position, we also separately fit the center positions 
of the VLBA and MERLIN data. Ultimately, we aligned 
the VLBA data to the MERLIN data assuming a common 
center for the R4 ring. This latter assumption, equiva¬ 
lent to an assumption of zero proper motion of the central 
protostar, does not affect the model per se, but it does 
provide an astrometric alignment between the data sets 
that can be checked for self-consistency. We will return to 
this astrometry check below. 

Fitting involved a non-linear least squares technique 
minimizing the % 2 difference between the model and data 
positions and velocities. Among the MERLIN data, many 
fainter maser spots clearly do not associate with the arcu¬ 
ate structure traced so clearly by the R4 masers. To re¬ 
move these unassociated masers objectively, we repeated 
the fitting procedure after clipping data for which the de- 
projected radius disagreed with the best-fit ring radius by 
more than 3 a , where cr was calculated by a quadratic 
sum of the model radius uncertainty and the data posi¬ 
tional uncertainty. The results of this modeling for all free 
parameters are listed in Tables 1 and 2, and the sky and 
kinematical projections are presented in Figures 6 and 7 
for the 1996 VLBA and 2000 MERLIN observations, re¬ 
spectively. The enclosed mass is 3.2 M Q ; a ZAMS star 
of this mass would be of spectral type A0 / B9 (de Jager 
& Nieuwenhuijzen 1987; Palla & Stahler 1993; Drilling & 
Landolt 2000). 

Table 1 

Results of the rotating ring model fit to the R4 

MASER DATA. 


Parameter 

Value 


Units 

RA offset 

-7.0 

± 

0.1 

mas 

Dec offset 

-27.4 

± 

0.2 

mas 

Inclination 

50 

± 

1 

degrees 

PA 

142 

± 

2 

degrees 

Sys. velocity 

-12.1 

± 

0.5 

km s -1 

Vout (1996) 

5.5 

± 

0.5 

km s -1 

fout (2000) 

5.3 

± 

0.5 

km s -1 

Central mass 

3.2 

± 

0.2 

M 0 


Table 1 Comments: Position offsets refer to the location of 
the center of the ring and are referenced to the position of 
the brightest maser: a(2000) = 22 h 56 m 1719807, 5(2000) 
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Major Axis Offset (mas) 



Fig. 6.— Results of the rotating ring model for the 1996 VLB A observations of (Torrelles, et al. 2001a). In each panel, the model is 
plotted as a solid line. The VLBA data are plotted as symbols, and the type of symbol depends on the epoch of observation as described in 
the legend of the topmost panel. Sky coordinates are referenced to the model ring center (see Table 1). Upper Left Panel : data and model 
projected on to the plane of the sky. The best fit sky offset to the ring center has been subtracted from the sky coordinates. Upper Right 
Panel : Radial velocity vs. offset along the ring major axis. Bottom panel : radial velocity vs. offset along the ring minor axis. Note that the 
minor axis offset coordinates have been deprojected for ring inclination. 
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Minor Axis Offset (mas) 


Fig. 7. — Results of the rotating ring model for the 2000 MERLIN observations. In each panel, the model is plotted as a solid line. Open 
circles mark data that were rejected as not belonging to the ring during the fitting process; solid circles mark data used in the fit. Sky 
coordinates are referenced to the model ring center (see Table 1). Upper Left Panel: data and model projected on to the plane of the sky. 
The best fit sky offset to the ring center has been subtracted from the sky coordinates. Upper Right Panel: Radial velocity vs. offset along 
the ring major axis. Bottom Panel: radial velocity vs. offset along the ring minor axis. Note that the minor axis offset coordinates have been 
deprojected for ring inclination. 
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= +62 o 01'49"429 (accurate to about 12 mas). The best 
fit radii are listed in Table 2. 

Table 2 

Best fit radii for the rotating ring model fits to 
THE R4 MASER DATA. 


Epoch Radius (mas) Radius (AU) 


1996.12 

21.6 

± 

0.2 

15.7 

± 

0.1 

1996.19 

22.2 

± 

0.2 

16.1 

± 

0.1 

1996.28 

23.3 

± 

0.2 

16.9 

± 

0.1 

2000.27 

38.7 

± 

0.1 

28.1 

± 

0.1 


Figures 6 and 7 show that the expanding, rotating ring 
models describe well the velocity gradient and spot distri¬ 
bution; certainly, this disk interpretation provides a better 
match to the maser spot kinematics than does the bow 
shock model. The positions of most of the maser spots 
agree to within a few mas of the model ring positions. The 
average proper motion of the R4 arc, based on the increas¬ 
ing radius of curvature between 1996 and 2000 (Table 2), 
is 3.9±0.2 mas yr -1 ( 13.3±0.1 km s” 1 ). The radial veloc¬ 
ities generally agree to within 5 km s -1 . The scatter of the 
velocity residuals might be attributed to substructure in 
the shockwave and turbulence in the post-shock gas; such 
scatter is observed in masers associated with outflows (e.g., 
Elitzur 1992). 

There remain, however, significant outliers in both po¬ 
sition and velocity. It is difficult to account for large po¬ 
sitional outliers in the context of this shockwave model. 
Masers lying within the ring radius for a given epoch may 
arise in clumps with longer post-shock relaxation times; 
masers outside the ring might be produced in a radiatively 
heated shock precursor region (Tarter & Welch 1986). Ei¬ 
ther explanation is unfortunately speculative and difficult 
to test owing to the small angular scale and high extinc¬ 
tion of the region. On the other hand, the outliers may 
arise from clumps not physically associated with R4-A. 
Cepheus A is an active star forming region with many ra¬ 
dio sources (Sargent 1977; Hughes & Wouterloot 1984b; 
Hughes, Cohen, & Garrington 1995); sight-lines to the R4 
region may well intercept more than one system of maser 
spots. 

We investigated whether the MERLIN positional out¬ 
liers might yet arise from the putative disk but at radii 
different from the shock ring. If the outliers are part of the 
disk, then their radial velocities should be compatible with 
the Keplerian rotation curve used to model the rings. We 
calculated the theoretical radial velocities at the positions 
of each of the R4 maser spots and subtracted the model ve¬ 
locities from the observed velocities. The resulting residual 
field is plotted in Figure 8. The velocity residuals of maser 
spots located on the rings are small, as would be expected 
from inspection of Figures 6 and 7. More interestingly, the 
velocity residuals of the maser spots not associated with 
rings are substantially reduced, a result arguing that most 
of the maser spots, whether on the ring for a given epoch 
or not, appear to be kinematically associated with R4. 

Two groups of R4 masers are glaring exceptions; these 
groups have been circled on Figure 8. Both groups have ve¬ 
locities that are too near systemic for their projected posi¬ 
tion onto the model disk, and, as such, show the largest ve¬ 
locity discrepancies with respect to the disk model. Their 


velocities are in better agreement with the range of ve¬ 
locities among the R4-B and R4-D masers, located south¬ 
east of the R4 rings and forming an apparent “tail” ex¬ 
tending away from the R4 arc (see Figure 3). Figure 9 
shows the kinematical distribution of the outlier maser 
spots with radial velocities in the range —15 km s _1 < 
Vlsr < —5 km s“ , i.e., the range of velocities spanned 
by the R4-B and R4-D masers. It is interesting to see that, 
with the ring masers removed, the outlier groups appear 
to align with a larger arc traced by the R4-B and R4-D 
masers and extending through the R4 ring. We can only 
speculate whether or not this arc is somehow associated 
with the R4 ring; it seems more likely that it is a separate 
maser arc that happens to lie near the same sight-line to 
the R4-A arc. Further observations may reveal whether 
these outlier groups, if they persist, participate in a sys¬ 
tematic proper motion with the R4-B and R4-D masers. 

4.3. Discussion of the Disk Model 

The models described above are necessarily simplifica¬ 
tions, but they provide the basic measurements of the ex¬ 
pansion of the R4 masers needed to deduce some prop¬ 
erties of the shockwave and the pre-shock gas. The disk 
model is somewhat limited by the assumption of spher¬ 
ical symmetry of the shockwave and cylindrical symme¬ 
try of the disk. Lacking additional data for the source of 
the shockwave, one might relax the condition of spherical 
symmetry and find a better fit to the data. For exam¬ 
ple, a more collimated outflow inclined into the plane of 
the disk might explain the lack of masers along the minor 
axis of the R4 maser distribution. The R4 maser spots 
sample the local kinematics unfortunately too sparsely to 
justify a more sophisticated model. It seems nevertheless 
likely that the shockwave must at least propagate directly 
along the projected major axis, and the properties of the 
spherical model apply to the shockwave along that axis. 

Referring to Table 2, the length of the semi-major axis 
expanded by 11 AU between 1996 and 2000, corresponding 
to an average ring expansion speed of 13 km s . The ex¬ 
pansion speed over the two months of VLBA observations 
is somewhat larger, roughly 30 - 40 km s^ 1 based on the 
proper motion between the 1996.12 and 1996.28 epochs, 
and so it seems the ring expansion has been decelerating. 
The cause of the deceleration may simply be mass-loading, 
analogous to the snowplow phase of a supernova shell, al¬ 
though scaled down in energy and momentum. 

There are unfortunately insufficent data to fit even a 
uniform deceleration model believably; the simplest mod¬ 
els (uniform deceleration and snowplow deceleration) have 
two free parameters, but there are only four data points. 
Nevertheless, to get an estimate of how much the ex¬ 
pansion speed may differ from the 1996 - 2000 average 
speed, we fit a snowplow model (after Dyson & Williams 
1997), appropriate for expansion into a dense medium, to 
the data in Table 2. We anchored the initial radius to 
15.7 AU (the best-ring radius at epoch 1996.12) and al¬ 
lowed the initial velocity to vary. The best-fit initial ve¬ 
locity is 41.3 ± 0.7 km s -1 , predicting an epoch 2000.27 
velocity of 7.2±0.7 km s -1 . These values should be viewed 
with caution, of course, because the present data do not 
sufficiently sample the ring expansion well enough to sup¬ 
port the snowplow model in particular. The snowplow 
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Fig. 8.— Residual analysis of the rotating ring model. Sky coordinates are referenced to the model ring center (see Table 1). The symbols 
marking the R4 maser positions are shaded according to their residual velocity ( vrsr — model in km s —1 ). Circles represent the MERLIN 
(2000) data, and squares mark the VLBA (1996) positions of the R4-A masers. The dotted line traces the best fit ring to the MERLIN data, 
and the smaller, gray ellipses surround the two spot groups with the largest velocity discrepancies with respect to the rotating ring model. 



150 100 50 0 -50 

RA Offset (mas) 


Fig. 9.— The MERLIN positions of the R4 masers within the velocity range 15 km s —1 < vlsr < 5 km s -1 . Maser spots whose kinematics 
are well described by the rotating ring model (velocity residuals < 5 km s — - 1 ) have been removed to emphasize the distribution of the outliers. 
The spots are shaded based on their radial velocity. The dotted line traces the best fit ring to the MERLIN data, and the darker ellipses 
surround the two spot groups with the largest velocity discrepancies with respect to the rotating ring model. 
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model does however show that the expansion speed at 
epoch 2000.27 may differ from the average by nearly a 
factor of two. 

Following the arguments of § 4.1, the ring expansion 
speed should be identically the shock speed. In the discus¬ 
sion to follow, we scale the shock speed to the average ring 
expansion speed, 13 km s , accepting that the instanta¬ 
neous shock speed may be somewhat lower in the later 
epochs. Otherwise, equating the ring expansion speed 
and the shock speed is valid provided the cooling time 
does not significantly vary with distance from the source 
of the shockwave. It has been proposed that the condi¬ 
tions necessary for H 2 O maser emission may be caused 
either by fast, dissociative shocks (Elitzur, et al. 1989) 
or slow, non-dissociative hydromagnetic shocks (Kaufman 
& Neufelcl 1996). We can rule out dissociative shocks for 
R4, which require vs > 50 km s~ . Our concern now is 
whether slower shocks through a protostellar disk could 
produce the observed high brightness temperatures. 

Kaufman & Neufelcl (1996) generated an array of shock 
models to calculate the efficiency of slow, C-type hydro- 
magnetic shocks. They defined the efficiency as the ratio of 
the maser luminosity output divided by the total mechani¬ 
cal energy provided by the shock: e sa t = 2L sat /A m noi-iVg, 
where L sat is the luminosity of the (saturated) water 
maser, A m is the surface area of the maser source, no is 
the pre-shock hydrogen density, and /i = 4.2 x 1CK 24 g is 
the mean mass per hydrogen atom. The maximum maser 
luminosity occurs where the product no£ sa t is a maximum. 
Using vs = 13 km s -1 and interpolating Figure 6 of Kauf¬ 
man & Neufeld, the maximum luminosity occurs for no = 
6.9 x 10' cm -3 , at which e sat = 4.0 x 10 -5 . The predicted 
luminosity is then L S at = 3.9xl0 25 (£/AU ) 2 erg s 1 where 
£ is the size of the maser cloud. The masers are probably 
beamed; we will make the simplifying assumption that the 
maser arises from a cylinder of cross-sectional diameter d. 
Taking into account the correction for cylindrical beam¬ 
ing, an observer in the path of the beam would infer an 
isotropic luminosity: 

L iso = 3.9 x 10 25 ^erg s ' 1 (4) 

VLBA measurements limit d < 0.4 AU (Torrelles, et al. 
2001 a), giving 

Liso > 2.4 x 10 26 er § S_1 ( 5 ) 

The flux density of the brightest R4 maser spot is 650 
Jy in a 0.2 km s -1 channel. Assuming isotropic emis¬ 
sion, the luminosity of any of the maser spots is therefore 
L 0 bs < 6.0 x 10 27 ergs _1 . The size of the maser cloud re¬ 
quired by the Kaufman & Neufeld (1996) model is found by 
equating the observed ( L 0 t , s ) and predicted ( Li SO ) maser 
luminosities: £ < 2.2 AU. It seems that the slow shocks 
model could plausibly produce the observed maser lumi¬ 
nosities insofar as the inferred path length is less than 7% 
the radius of the disk. For comparison, the gas scale height 
derived assuming hydrostatic equilibrium (and neglect¬ 
ing self-gravity) is H ~ c s R 3/2 / ( GM) 1/ 2 ~ 5.4T 3 1/2 AU, 
where T 3 = T/1000 K (Frank, King, & Raine 1992). Opti¬ 
cal and infrared imaging of circumstellar disks in a variety 


of environments measure scale heights ranging from a few 
AU to tens of AU at comparable distances from the cen¬ 
tral star (e.g., Beckwith & Birk 1995; Stapelfeldt et al. 
1998; Heap et al. 2000; Stapelfeldt 2000; Stapelfeldt et al. 
2000 ). 

The parameters of the shock model, taking into ac¬ 
count the properties inferred from the Kaufman & Neufeld 
(1996) model, are sufficient to allow us to estimate the 
mass of the circumstellar disk. Assuming constant scale 
height, uniform density, and cylindrical symmetry, the disk 
mass is: 

Mdisk ~ 0.97 x 10 7 cm _ 3 ) ( 2 8 Au) (^2Au) M ® 

( 6 ) 

where we have normalized R to the radius of curvature 
based on the MERLIN observations (Table 2). Note that 
the pre-shock density is optimized for the most luminous 
maser emission, as discussed above. Such a low disk mass 
is reminiscent of Vega-like debris disks, whose masses are 
in the few range (e.g., Sylvester & Skinner 1996). The 
disks surrounding young stars and protostars are measured 
to be ~ 0.01 a few Mq (e.g., Hillenbrand et al. 1992; 
Chandler & Richer 1999 & references therein). Whereas 
it is tempting to infer that the R4 system might be more 
evolved than the surrounding protostars of the Ceph A re¬ 
gion, it seems just as likely that the gas associated with 
R4 is very clumpy. Masers might then arise from relatively 
low density gas but are quenched in higher density regions. 
Nevertheless, the mass estimate is at least self-consistent 
with our original assumption of negligible disk mass. As 
the ring expands and traces out the rotation curve, future 
epochs of maser observations may place better dynamical 
constraints on the disk mass. 

Based on the kinematic model fitting, the expansion ve¬ 
locity of the R4 maser spots is roughly 5 km s ” 1 (Table 1), 
which we interpret as the post-shock velocity. The differ¬ 
ence between the shock velocity and post-shock velocity is 
somewhat surprising: a strong shock should produce post¬ 
shock velocities 3ns/4 < v' < vs (McKee & Hollenbacli 
1980; Dyson & Williams 1997). Even allowing for decel¬ 
eration, the velocity difference remains an issue for the 
1996 epochs, during which the expansion speed is of order 
40 km s _1 , eight times the post-shock speed. Weakening 
the shock to M ~ 1.5, referenced to the average shock ve¬ 
locity, would drop the post-shock velocity to 5 km s _1 but 
implies pre-shock temperatures T > 10 4 K in order to raise 
the sound speed sufficiently. The pre-shock gas should be 
molecular to produce post-shock masers (Elitzur, et al. 
1989; Kaufman & Neufeld 1996), and so we can rule this 
option out. 

Another explanation of the low post-shock velocity is 
that the magnetic pressure of the post-shock gas may dom¬ 
inate the gas pressure (see Liljestrom & Gwinn 2000 for 
an identical analysis of W49N). Magnetic field pressure 
inhibits compression of the post-shock gas, which, by con¬ 
servation of momentum, increases the speed of the gas 
leaving the shock. To the observer, then, the effect is to 
reduce the post-shock velocity. In the limit where mag¬ 
netic field pressure dominates gas pressure, the magnetic 
pressure balances the ram pressure of the pre-shock gas 
(Elitzur, et al. 1989; Kaufman & Neufeld 1996). The so- 
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lution for the post-shock magnetic field strength in terms 
of the pre-shock gas density 717 = n/ 10 7 cm -3 and shock 
velocity Vs is: 


B' = 31 


vs 


13 km s 


-1 


0.5 


mG 


( 7 ) 


In the frame of the shock, the pre-shock magnetic field is 
given by (Hollenbach & McKee 1979; Liljestrom & Gwinn 
2000 ): ‘ 


n° 7 5 mG, (8) 

(9) 

where parallel-line subscripts indicate components parallel 
to the shock front. A magnetic field strength of ~ 30 mG 
is plausible in comparison with more direct measurements 
of magnetic field strengths in similar astrophysical envi¬ 
ronments. For example, the magnetic field strengths as¬ 
sociated with the OH masers of Cepheus A are a few mG 
(Wouterloot, Habing, & Herman 1980; Cohen, Brebner, 
& Potter 1990). Magnetic field strengths of 20 - 100 mG 
have been measured by observations of Zeeman splitting 
of H2O masers found in other star-forming regions (Fiebig 
& Guesten 1989; Sarma, Troland, & Romney 2001). Mea¬ 
suring Zeeman splitting of the H30 a recombination line, 
Thurn & Morris (1999) measured a field strength of 22 mG 
in the corona of the circumstellar disk of MWC 349. 



5. PROPER MOTIONS OF THE R2, R3, AND R5 MASER 
REGIONS 

The four-year separation between the VLBA and MER¬ 
LIN observations affords an opportunity to study proper 
motions over the entire HW2 region. Unfortunately, the 
1996 VLBA observations and the present, 2000 MERLIN 
observations are referenced to different maser spots, and 
furthermore the 1996 VLBA data are not phase-referenced 
to a fixed calibrator position. We can use, however, the 
proper motions of the R4 masers to calculate a boot¬ 
strapped astrometry between the data sets. Based on the 
increasing radius of curvature of the R4-A masers (and ir¬ 
respective of the actual nature of the R4 masers, whether 
disk or outflow or other), it seems clear that the R4-A 
region is expanding away from some common center. 

As a first guess, we assume that proper motion of the 
R4 expansion center is negligible. Figures 10-12 plot the 
relative positions of the MERLIN and VLBA data for the 
regions R1-R5. The MERLIN positions of the R2 and R3 
maser spots are displaced by roughly 5 mas to the west of 
the positions predicted by the proper motions measured 
by Torrelles et al. (2001a), but the declination alignment 
is somewhat better. From inspection of Figures 4 and 8 of 
Torrelles et al., the proper motion uncertainties are prob¬ 
ably of order 1 mas yr , which propagates to ~ 4 mas 
between the MERLIN and VLBA data sets. It is there¬ 
fore unclear whether the displacement might owe to proper 
motion of the R2 and R3 spots, proper motion of the R4 
spots, which were used as the astrometric reference, or 


both. Follow-up, phase-referenced MERLIN observations 
should answer this question. 

Torrelles et al. (2001b) emphasized the circular symme¬ 
try of the R5 region and its apparent radial expansion over 
the 2 months spanned by the VLBA observations. We re¬ 
measured the proper motion of the R5 structure by fitting 
circles separately to each VLBA epoch. The R5 ring ex¬ 
pands at a rate of 2.5 ± 0.1 mas yr -1 (8.6 ± 0.3 km s -1 ), 
and the expansion center moves at 1.4 ± 0.1 mas yr -1 (4.8 
± 0.3 km s _1 ) towards P.A. 126° ± 6°. For comparison, 
Torrelles et al. reported an expansion speed of ~9 km s _1 
and center motion of ~ 6 km s _1 into PA ~ 143°. 

We calculated the location and geometry of the R5 arc 
as it should have appeared for the (2000) MERLIN obser¬ 
vations based on the expansion measured by the (1996) 
VLBA observations; the results are plotted in Figure 12. 
Six out of the ten R5 maser spots detected by MERLIN 
fall within the predicted range of sky offsets, and four are 
displaced south of the prediction band by ~ 3 mas. We 
judge this agreement to be very good, considering the un¬ 
certainty in predicting four years of proper motion on the 
basis of two months of monitoring. Based on the accept¬ 
able agreement for the R2, R3, and R5 masers, the original 
assumption for the astrometry bootstrap, namely, negligi¬ 
ble proper motion of the R4-A curvature center, seems 
reasonable. 

The argument can be reversed to constrain the proper 
motion of the R4 masers. Based on the ~ 5 mas western 
displacement of the R2 and R3 masers, a < 1.2 mas yr -1 . 
The southern displacement of the R5 masers gives 5 < 
0.8 mas yr -1 . Taking these values as a limit on the mo¬ 
tion of the astrometric reference, the proper motion of the 
R4 expansion center must be less than ~ 1.4 mas yr -1 (5 
km s -1 ). 


6. SUMMARY AND DISCUSSION 

The simplest explanation for the expansion of the R4 
maser arcs appears to be a “slow,” C-type liydromagnetic 
shockwave propagating through a rotating, circumstellar 
disk surrounding a forming or young A0 / B9 star. Owing, 
in part, to the high extinction towards this region, the cen¬ 
tral star has yet to be identified at any waveband (HW2 
itself has only been identified in radio continuum). It is 
therefore difficult to measure or constrain any other prop¬ 
erties of the young star and its associated shock except for 
what we have learned from the maser kinematics. 

The YSO IRAS 21391+5802 shows a similar, but smaller 
(~ 2 AU diameter) ring of masers, which presumably sur¬ 
round the young star (Patel et al. 2000). In that source, 
the masers seem to occur at the dust condensation ra¬ 
dius of a radial outflow from the YSO. As the molecu¬ 
lar gas flows outward from the dust condensation front, 
the gas cools, and the masers associated with that mate¬ 
rial fade out. The maser ring persists as fresh material 
enters the dust condensation front. This model predicts 
that the diameter of the maser spot distribution should 
not change, however, in contrast to the appearance of the 
R4 maser arcs, and furthermore the model predicts no 
strong velocity gradient around the ring. The R4 maser 
arcs are furthermore non-circular and much larger, extend¬ 
ing well outside the dust condensation radius (the masses 
of the two YSOs are similar). It seems that, in contrast 
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Fig. 10.— Comparison of the MERLIN and VLBA data for maser regions Rl, R2, and R3. MERLIN detected no maser emission that 
could be clearly associated with Rl. The gray bands trace the predicted 2000 positions of the Rl and R3 masers based on the proper motion 
analysis of Torrelles et al. (2001a). The alignment between the MERLIN and VLBA data sets is based on the assumption that the proper 
motion of the expansion center of the R4-A masers is negligible. 



Fig. 11.— Comparison of the MERLIN and VLBA data for the R4 maser regions. The alignment between data sets is based on the 
assumption that the proper motion of the expansion center of the R4-A masers is negligible. 
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Fig. 12.— Comparison of the MERLIN and VLBA data for the R5 maser region. The alignment between data sets is based on the 
assumption that the proper motion of the expansion center of the R4-A masers is negligible. The VLBA data are plotted as open symbols, 
and the MERLIN data are plotted as filled circles. The curved lines trace the region in which the MERLIN data should appear based on 
the 6.2 ± 1.8 mas yr —1 proper motion measured from the three closely-spaced VLBA epochs. The solid line traces the nominal, predicted 
location, and the dotted lines mark the 90% confidence prediction bands. 


to IRAS 21391+5802, the R4 arcs cannot be explained as 
a standing feature of a steady outflow; the maser proper¬ 
ties are better described by gas responding to a shockwave 
radially propagating through a rotating medium. 

The origin of the shockwave is unclear. One possibility 
might be colliding spiral shocks, or other arcuate struc¬ 
ture, forming as the result of disk instabilities; for example, 
Durisen et al. (2001) proposed this scenario to explain the 
occurrence of methanol masers in protostellar disks. Disk 
and pattern rotation would alter the shape of the maser 
arc over time, but it does not seem likely that spiral shocks 
would produce the observed, outward propagation of the 
R4 maser arcs, nor should they produce the increasing ra¬ 
dius of curvature. Spiral shock models also require that 
the disk is viewed nearly edge-on for significant amplifica¬ 
tion (Maoz & McKee 1998), but it seems that the R4 arcs 
are viewed at an intermediate inclination. 

A more promising explanation is provided by the work of 
Tscharnuter, Boss, & collaborators (Morfill, Tscharnuter, 
& Volk 1985; Tscharnuter 1987; Boss 1989). Using hy- 
drodynamical simulations, they have demonstrated that a 
collapsing cloud core can become unstable to large oscil¬ 
lations owing to the thermodynamics of molecular hydro¬ 
gen dissociation and reassociation. As the protostar core 
forms, this “hiccup” instability can drive AU-scale, radial 
outflows of order 10 km s -1 (Boss 1989), comparable to 
the outward proper motion of the R4 maser arcs. The hic¬ 
cups recur as infalling material drives the protostar again 
to instability. Balluch (1988) argued that this cycle of in¬ 
stability proceeds for at least hundreds of years (i.e., at 
least as long as the simulated duration of the numerical 


models). We speculate that perhaps the disk associated 
with the R4 masers may be responding to this sort of in¬ 
stability of the protostar. It is difficult to evaluate the en¬ 
ergetics involved based on the simulations that have been 
published to date, and we are unaware of any work model¬ 
ing the impact of the hiccup instability on a surrounding 
protostellar disk; such analysis is beyond the scope of the 
present work. In addition, the hiccup models require fairly 
high infall rates, which we do not have the data to justify 
in the case of the R4-A masers. As was pointed out by 
Torrelles et al. (2001a; 2001b) in their discussion of the 
neighboring maser source R5, the origin of disk masers is 
yet poorly understood. Our shock interpretation argues 
for more work towards understanding instabilities in pro¬ 
tostars and the impact of such instabilities on circumstellar 
disks. 
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